fx=[0 1 1 1 0];
subplot(3,2,1),imshow(fx),xlabel('(b) 原始图像');
subplot(3,2,2),stem(fx,'.'),axis([0 6 0 2]),xlabel('(c) 原始图像二维示意图');

Fu=fft2(fx);
ABSFu=abs(Fu);
subplot(3,2,3),imshow(ABSFu,[]),xlabel('(f) 傅里叶幅度谱图像');
subplot(3,2,4),stem(ABSFu,'.'),axis([0 6 0 4]),xlabel('(g) 傅里叶幅度谱二维示意图');

ShiftABSFu=fftshift(ABSFu);
subplot(3,2,5),imshow(ShiftABSFu,[]),xlabel('(i) 对称傅里叶幅度谱图像');
subplot(3,2,6),stem(ShiftABSFu,'.'),axis([0 6 0 4]),xlabel('(j) 对称傅里叶幅度谱二维示意图');